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Abstract. In this paper basic differential invariants of generic hy- 
7~-j ' perbolic Monge-Ampere equations with respect to contact transforma- 

^Z. ', tions are constructed and the equivalence problem for these equations 

^T; ■ is solved. 

U 

< 

Q\ . 1. Introduction 

With this paper we start a systematic study of differential invariants of 
Monge-Ampere equations, with the objective of the classification problem, 
^y^ I methods of integration and other applications. Complete proofs of the results 

pj ■ announced in [14] are presented. We are interested in the classical case of 

two independent variables. The Monge-Ampere equations equations merit 
a special attention due to a large spectrum of various applications, first of 
all, in differential geometry and mathematical physics. Moreover, they form 
a natural testing area for new methods emerging in the modern theory of 
nonlinear PDE's. 



> 

00 

^r^ . In spite of more than 200 years of history of Monge-Ampere equations 

^D I and numerous publications devoted to them it would be an exaggeration to 

say that their nature is well understood. An important success was establish- 

\^ , ing the existence and uniqueness theorems by Lewy and others (see [10, 3] 

^^ ' for local aspects and [20] for global ones). The classical Monge integration 

method was modernized by Matsuda [15, 16] and Morimoto [17], etc. Our 

interest in differential invariants is motivated not only by the classification 

fi ■ problem but, no less, by hopes that they could illuminate many aspects of 

^ . the theory of Monge-Ampere equations. 

According to [22] (see also [1]) scalar differential invariants provide a key 
}_j ■ to solving the classification problem for any kind of geometrical structures. In 

5t ! fact, geometrical structures of a given type are classified by solutions of a nat- 

urally associated classifying (differential) equation, which describes "family 
ties" connecting the corresponding scalar differential invariants. More ex- 
actly, scalar differential invariants are smooth functions on the classifying 
diffiety, which is the infinite prolongation of the classifying equation. This 
diffiety has, generally, singularities and its singular strata classify those geo- 
metrical structures that possess nontrivial symmetries. Each of these strata 
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is also an infinitely prolonged differential equation in a lesser number of inde- 
pendent variables. For instance, homogeneous structures correspond to the 
zero-dimensional case. So, the classification problem consists of a complete 
description of all strata composing the classifying diffiety and, therefore, 
involves a complete symmetry analysis of the geometric structures under 
consideration. The interested reader will find an illustration of the above 
said in [23] where plane 3-webs, a rather simple geometrical structure, is 
considered. 

The classification problem for Monge-Ampere equations dates back to 
Sophus Lie. For modern proofs of Lie's theorems, classification problems for 
various strata of Monge-Ampere equations see, e.g., [13, 6, 7, 8, 9, 18] and 
references therein. In this paper we interpret a hyperbolic Monge-Ampere 
equation as a pair of 2-dimensional, skew-othogonal non-lagrangian subdis- 
tributions of the contact distribution on 5-dimensional contact manifolds. 
Another approach to these equations was developed by V.V. Lychagin in 
[11, 12]. We look for differential invariants of Monge-Ampere equations, not 
only scalar, with respect to the group of contact transformations. Here we 
limit ourself to the case of generic hyperbolic equations. This is motivated 
by two reasons. First, the study of singular strata benefits much from the 
knowledge of the generic one. Second, for the hyperbolic equations differen- 
tial invariants are easier visible due to the existence of bicharacteristics. 

Differential invariants found in this paper give a solution of the classi- 
fication problem for generic hyperbolic equations. This solution requires a 
substantial computer support in analysis of concrete cases and a further 
work is necessary to improve its efficiency. 

Differential invariants for elliptic and parabolic Monge-Ampere equations 
can be obtained more or less straightforwardly by following the approach 
developed in this paper. This and a the study of singular strata will be the 
subject of subsequent publications. 

2. Preliminaries 

Below, all manifolds and maps are supposed to be smooth. By [f]p , k = 
0, 1, 2, ... , cx), we denote the /c-jet of a map / at a point p. M stands for the 
field of real numbers, and M" for the n-dimensional arithmetic space. 

2.1. Jet bundles. Here we recall necessary definitions and facts about jet 
bundles, see [4, 5]. 

Let M be an n-dimensional manifold, E an n + m-dimensional manifold 
and 

TT : E — > M . 

a fiber bundle. By 

TTk:j''7T^M, TTi, : [S]'; ^ p , fc = 0,1,2,... 

we denote the bundle of all /c-jets of sections of n. For any I > m > 0, the 
natural projection is defined as 

'^1,'m ■■ J IT ^ J^vr , TTi^jn ■ [S]p ^ [5]™ . 

Any section S" of vr generates the section jkS of the bundle tt^ by the formula 

ifc 



j,S:p^[Sfp 
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Put 

Let 6'fc_|_i be an arbitrary point of j'^'+^vr, 9k = 7rfc+i,fc(^fc+i), and T0^{J^it) 
the tangent space to J^t: at the point 6^. Then 6*^+1 defines the subspace 
Kq^^-^ C T0^{J^tt) by the formula 

Clearly, O^^i is identified with Kg . It is easy to prove that 

r,,(jV)=Ke,+i©T,,(7r,-i(p)). (1) 

Consider all submanifolds of the form Lg containing 9k- Subspace spanned 
by their tangent spaces T0^,{Lg) is denoted by Q{9k) and it is called the 
Cartan plane at 9k. The distribution 

Qk '■ Ok ^^ Q{9k) 

is called the Cartan distribution on J^tt. The distribution Sfc, /c > 1, can be 
defined as the kernel of the Cartan form 

Uk = pr2 O (TTk^k-l)* , 

where pr2 : Tg^_^(J^^^7r) -^ Tg^_^{-K^_^{p)) is the projection generated by 
direct sum decomposition (1). 

2.2. The contact structure. Consider the trivial bundle 

r:M^xM — > M^ t : {x,y,z) ^ {x,y) . 

By X, y,z, p = Zx, q = Zy, r = z^x, s = Zxy, t = Zyy we denote the standard 
coordinates in J'^t. 

The Cartan distribution Ci on J^t is identical to the contact structure 
on J^T. The corresponding contact 1-form Ui has the canonical form 

Ui = dz — p dx — q dy . 

in the standard coordinates. 

A diffeomorphism ip : J^t ^ J^t is called a contact transformation if 
it preserves the Cartan distribution. Obviously, a diffeomorphism (^ is a 
contact transformation iff there exist a nowhere vanishing function A such 
that 

cp*{Ui) = XUi. 
Any contact transformation ip can be lifted to the diffeomorphism 

^(1) : jV -^ J\ 
by the formula 

99« : 92 = Kg, ^ MKe,) = 02 = (/^l')(e2) • 

If if is defined on an open set V C J^t, then (pr is defined on an open, 
everywhere dense subset of Tq^-^^ {V). 

A vector field Z in J^t is a contact vector field if its flow ipt consists of 
contact transformations. Clearly, Z is a. contact vector field iff there exist a 
function A such that 

Lz{Ui) = \Ui, 
where Lz is the Lie derivative with respect to Z. 
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There exists a natural one-to-one correspondence between the set of ah 
contact vector fields in J^t and the set of all functions in J^t. It is defined 
by the formula 

Z ^ f = Z JUi . 

The function / = Z JUi is called the generating function of the contact 
vector field Z. The contact vector field Z corresponding to / is denoted by 
Zf. In standard coordinates, the field Zf is given by the formula 

Zf = -/p^-/.^ + (/-^/p-'?/.)^ + (/-+P^)^ + (/. + 9/.)|- (2) 

2.3. Operations over vector-valued forms. Let M be a smooth n-di- 
mensional manifold, A*(M) the C°°(M)-module of i-forms on M and D{M) 
the C°°(M)-module of vector fields on M. Let a G A^{M), (3 G A''(M), 
and X, y G D[M). Then the Frolicher-Nijenhuis bracket |- , •] of the vector- 
valued forms a® X and /3 y is defined by the formula 

la®X,[3®Y} 

= a A /3 [X, y] + a A X(/3) ®Y - Y{a) A P X 

+ {-l)''da A (X J/3) y - (-l)^(y Ja) A dp X , 

see [2]. 

The contraction J of forms a fS> X and (3 fS>Y is defined by the formula 

(a X) J (/3 y) = Q A (X J/3) y . 

2.4. Projectors and their curvatures. The following simple construction 
allows one to associate a vector valued 2-form with a projector. Namely, let 
P,Q £ D{M) be endomorphisms of the C°°(M)-module D{M) such that 
QP = 0. Then 

nQXX,Y)=Q[P{X),P{Y)], X,YeD{M), (3) 

obviously, is skew-symmetric and C°°(M)-bilinear, i.e., a vector valued form. 
More precisely, it takes values in Im Q C D{M). If P : D{M) -^ D{M) is a 
projector, i.e., P^ = P, then the associated curvature form of P is defined 
to be 

%p = ni_p,p (4) 

with I = id£,(Af). 

3. Hyperbolic Monge-Ampere equations 

3.1. Monge— Ampere equations. The Monge-Ampere equation is a par- 
tial differential equation of the form 

A \ZxxZyy — Z^y) + AZxx + DZxy + ^ Zyy + iJ ^ (J , (5j 

where x, y are independent variables, z is a dependent variable, Zxx = 
d'^z/dx^, Zxy = d'^z/dxdy, Zyy = d'^z/dy'^, and coefficients N, A, B, C, 
D are functions of x, y, z, Zx = dz/dx and Zy = dz/dy. 

We identify equation (5) with the submanifold £ of the jet bundle J'^t 
determined by the equation 

N{rt-s'^) + Ar + Bs + Ct + D = 0. (6) 
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Obviously, 

T2,i(£)=JV. 

Let 02 G £, T2,i(02) = Gi-, and Fq^ be the fiber of tlie projection r2,i over 
the point 9i € J'^t. Then the subspace 

Smbl02 £ = Te2£ n Te^Fe^, 

where T02£ is the tangent space to £ at O2 is called the symbol of the equation 
£ at the point ^2 G £• In terms of standard coordinates, Smblgj £ is described 
by the linear equation 

N{tf + ri- 2ss) + Af + Bs + Ci=f), (7) 

where f, s, t are the standard coordinates in T0,^F0^ generated by the stan- 
dard coordinates on J'^t. 

A point 02 G £■ can be elliptic, parabolic, or hyperbolic. To introduce 
these notions, let us consider a one-dimensional subspace P C S(^i) such 
that {ti)^P / 0. By definition, put 

l{P) = {02eFe,\PcKe,}- 
The submanifold 1{P) is called a 1-ray. In terms of standard coordinates, let 
01 = {x,y,z,p,q), P = {v) and 

Odd d d 

dx dy dz dp dq 

Then {ti)^P 7^ means that 

(Ci,C2)/(0,0), (9) 

V £ C(^i) means that 

f^ = CiP + C2q, (10) 

and P C Kg^ means that 

r?i = Cir- + C2S , 

V2 = Cis + C2t, ^ ^ 

where r, s, t are the standard coordinates of 02 in the fiber Fg^. From system 
(11), we see that 1{P) is an affine straight line in Fg^^. By i.Q^{P) we denote 
the tangent space Tg^l{P) to 1{P) at the point 62 G K-f)- We call it a i-ray 
subspace. In terms of the standard coordinates r, s, t in Tg^Fg-^, vectors of 
^02 (P) satisfy 

Ulf + C2S = 0, 

|ci5 + C2t = o, 

Obviously, £^2 i^) is spanned by the vector 

{f,s,i) = {cl-CiC2,c!)- (13) 

Taking into account (9), we observe that all 1-ray subspaces form the cone 

Vg,={fi-S^=0} 

in the tangent space Tg^ Fg^ . This cone is called the cone of singular square 
forms. Obviously, the intersection Smble2 £ H ^512 is either zero, or a single 
1-ray subspace, or two 1-ray subspaces. Correspondingly, the point 02 G 8, 
is then called elliptic, parabolic or hyperbolic. It is not difficult to prove that 
a contact transformation takes an elliptic, parabolic, or hyperbolic point to 
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an elliptic, parabolic, or hyperbolic point, respectively. The equation £ is 
called elliptic, parabolic or hyperbolic if all its points are elliptic, parabolic 
or hyperbolic, respectively. In this work, we consider hyperbolic Monge- 
Ampere equations only. It is easy to see that £ is hyperbolic iff its coefficients 
satisfy the condition 

A = B^ -4AC + 4ND>0. (14) 

3.2. Skew-orthogonal distributions. Following [21], we show that a hy- 
perbolic Monge-Ampere equation is equivalent to a pair of skew-orthogonal 
two-dimensional distributions in the Cartan distribution on J^t. 

Let 6i be an arbitrary point of J^t. By Qgj we denote the union of all 
one-dimensional subspaces P of C(^i) such that t*P 7^ and the 1-ray 1{P) 
is tangent to £ at least at one point. 

Proposition 3.1. Let 8, be a hyperbolic Monge-Ampere equation. Then Qg^ 
is the union of two-dimensional subspaces ^'^(^i) and 'D'^{9i) of the Cartan 
plane G{6i), so that 

(1) e(ei)=!Di(0i)©D|(0i), 

(2) ©^(^i) and 1)'^{6i) are skew-orthogonal with respect to the symplectic 
form dUi = dx f\dp + dy f\dq on C. 

Proof. We prove this proposition for Monge-Ampere equations such that 
N ^ 0. The proof for A^ = follows from the fact that every Monge- 
Ampere equation can be transformed to one with A^ 7^ by an appropriate 
contact transformation. 

Let V £ Qg-^ and P = {v). The condition for 1{P) to be tangent to £ can 
be written in the following way. We can assume that v is of the form (8). 
Then the vector of fiber coordinates (Cl) ~CiC2) Ci) is tangent to 1{P). Now 
using (7) we deduce that 1{P) is tangent to £ iff 

N{rCf + 2<iC2 + tCl) + ACl - BC1C2 + CC? = . 

Taking into account that the coordinates Q and rji of v are connected by 
equations (11), we reduce this equation to the form 

NiCim + C2m) + Acl - i?CiC2 + ccl = o . (is) 

Then in view of (9) we assume that Ci 7^ (the case C2 7^ is analogous). 
Then from (11) we get 

r = -p2{riiC.i - ^2^2 + Ch) , s = —{r]2 - C,2t) . 

Substituting these expressions for r and s in equation (6) and taking into 
account equation (15), we obtain the equation 

Nr,l + {AC2 - BCi)m - ACm - DCl = 0. (16) 

Solving the system of equations (15) and (16) with respect to r/i and 172, 
we find 

^ (^ T VA)C2 - 2CCi ^ {B± VA)Ci - 2AC2 

VI 2iV ' "^^ 2N 



MONGE-AMPERE EQUATIONS 

Finally, in view of (10), we see that 

f d d C d B±VA d 



dx dz N dp 2N dq 



^ . d d BtVA d Ad. ,^^, 



This proves that Qg^ = (Xi,X2) U {X^,X4) with 

X-— —-—— B-^ d 

dx dz N dp 2N dq ' 

_a d B + ^ d Ad 

dy dz 2N dp N dq' 

X-— —-—— -^ + ^ ^ 

dx dz N dp 2N dq ' 

a d B-^ d Ad 

dy dz 2N dp N dq 



(18) 



Put 



v\{e^) = {Xi,X2), vUei) = {Xs,x,) 



Now it is straightforward to verify that subspaces ©^(^i) and D|(^i) are 
skew-orthogonal and ©^(^i) n 'D'^{9i) = {0}. This completes the proof. D 

From (18) we see that for a Monge- Ampere equation such that N ^ 0, 
the map ri* projects ©^(^i) and D|(^i) onto the tangent space to the base 
of the bundle r without degeneration. 

It should be noted that if A^ = (that is, if £ is a quasilinear second 
order PDF), then the projections Ti^(T>j^{6i)^ and ri*(l'|(^i)) are one- 
dimensional. 

Thus an arbitrary hyperbolic Monge- Ampere equation generates two 2- 
dimensional skew-orthogonal subdistributions of the Cartan distribution Ci 
in J^T. 

Proposition 3.2. Let £. be a hyperbolic Monge-Ampere equation. Then 
92 ^ ^ if and only if one of the following equivalent conditions holds: 

(1) iiTgg r\'D\^{9i) is 1- dimensional, 

(2) Kq,^ n 1)1(^1) is 1- dimensional. 

Proof. As in the proof of Proposition 3.1 one can assume that A^ 7^ 0. 

Let 6*2 G £. Then Smblg^ £ n V^^ = ^e2((^')) U ^e^iiv)), where ioiiiv)) and 
i02{{v)) are different straight lines and, so, vectors v and v are independent. 
They are skew-orthogonal, since K^^ is a Lagrangian plane in S(^i) and, by 
definition of Q^^, v, v £ Qg^. This means that Kg.^ intersects planes ©^(^i) 
and I'|(6'i) along (v) and {v), respectively. 

Let 6*2 be a point of J'^t such that Kg^, intersects the plane ©^(^i) along a 
straight line, that is, Kg.^ riT)^^{9i) = {v). By substituting coordinates ??i,r/2 
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of the vector v given by formula (17) into eq. (11), we obtain 



IN P \ N P 

By hypothesis this system is of rank 1 (cf. (9)) and hence its determinant 
is zero. Now it remains to note that this is exactly equation (6) and, so, 
^2 G £• The case of T)\{d\) differs only by the sign at ^A. D 

An important consequence of this proposition is that a hyperbolic Monge- 
Ampere equation £ is completely determined by one of the associated dis- 
tributions Dg, i = 1, 2. 

Thus, every hyperbolic Monge- Ampere equation £ is naturally equivalent 
to a pair of 2-dimensional, skew-othogonal non-lagrangian subdistributions 
Dg, D| of the Cartan distribution Si in J^t. In particular, the equivalence 
problem for hyperbolic Monge-Ampere equations with respect to contact 
transformations may be interpreted as the equivalence problem for pairs of 
2-dimensional, skew-orthogonal non-lagrangian subdistributions of Ci with 
respect to contact transformations. 

3.3. Bundles of Monge Ampere equations. From now on we put M = 

3.3.1. Bundles of hyperbolic Mange-Ampere equations. Let £ be a Monge- 
Ampere equation (5). It is identified with the section 

5£ : X ^ [N{x) : A{x) : B{x) : C{x) : D{x)] 

of the trivial bundle 

p : MP^ X M — > M , ([p° : p^ : p2 :/:/], x) H^ X , 

where MP^ is the 4-dimensional projective space. Obviously, this identifica- 
tion is a bijection of the set of all Monge-Ampere equations onto the set of 
all sections of p. 

Consider the open subset E of the total space of p defined by the condition 
(14), i.e.. 

Clearly, the section S'g corresponding to a hyperbolic Monge-Ampere equa- 
tion £ takes values in E. Thus we can define the bundle of hyperbolic Monge- 
Ampere equations by the formula 

7r = p\j^:E — >M, {[p^ : p^ : p^ : p^ : p\ x) ^ x . (19) 

We use local coordinates x, y,z,p,q, n^, . . . , u^ in the total space E of tt, 
where x,y,z,p,q are the standard coordinates on M, while the coordinates 
u^, . . . , u'^ on the fibres of tt are defined as follows. Consider the affine 
hyperplane in M^ defined by the equation p^ = 1. It generates the local 
chart in E 

[1 : / : p2 : / : p^] ^ {p^ , p^ , p^ , p'^) ■ 
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Following formulas (18), we introduce the local coordinates u^, . . . , u^ along 
the fibres of tt by 

1 3 2 P^ - ^ 3 P^ + "^ 4 1 /onN 

« = -r , u = '-—^ , u^ = '-—^ , u^ = -p^ , (20) 

where A = (p^)^ — Ap^p^ + Ap'^. 

These coordinates extend to the standard coordinates x,y,z,p,q,u^,ul., 
Uy, ul,Up, Uq, . . . , «Jr, . . . , on j'^TT, used in this paper until we replace them 
with a more convenient set in Sect. 4.3. 

3.3.2. The lifting of contact transformations. Let (p he a contact transfor- 
mation defined in M. Then (p transforms any Monge- Ampere equation £ to 
another Monge- Ampere equation £. In other words, <p induces a transfor- 
mation of the corresponding sections S'g i— > Sp and, consequently, a diffeo- 
morphism ip^^> of the total space of n such that the diagram 

E -^ — > E 



M > M 

is commutative (in the domain of ^^^'). The diffeomorphism ip^^' is called 
the lifting of 99 to the bundle tt. 

The diffeomorphism ip^^\ in its turn, can be lifted to a diffeomorphism 
ip^ ' of J TT by the formula 



^('H[st) = [^(^Kso^~X^^^. 



Obviously, for any I > m, the diagram 



jl^ _^ jl^ 



J'^TT > J'^TT 



is commutative (in the domains of 99' -*). The diffeomorphism Lp^ > is called 
the lifting of pi to the jet bundle J tt. 

3.3.3. The lifting of contact vector fields. Let Z he a contact vector field in 
M and let (pt he its flow. Then ipl defines a vector field Z^^' in J'^tt. This 
field is called the lifting of Z to J^tt. Obviously, 

(7r;,^),(Z«)=Z(™), oo>/>m>-l, 

where Z^^^' = Z. 

It is not difficult to see that the map 

Z > — > Z^^'^ 

is a homomorphism of the Lie algebra of all contact vector fields onto the 
Lie algebra generated by all vector fields of the form Z^ ' . 

The local expression of Z^^' can be found as follows. First, change the 
notation by putting x"^ = x, x'^ = y, x^ = z, x'^ = p, x^ = q. Recall that 
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the operator Dj of total derivative with respect to x^ in Joq is given by the 
formula 

|o-|>0 i=l '^ 

The operator of evolution differentiation corresponding to a generating func- 
tion ij{Z) = {iIj^{Z), . . . ,^p'^{Z)Y is defined by the formula 



\a\>Oi=l '' 

where a = {ji . . . >} , ^o- = Dj^ o . . . o Dj^ and i'iZ) is defined as follows. 
Let S* be a section of vr defined in the domain of Z, 6i = [S]\, and 
X = 7ri(0i); then 

^l,{z){e,) = ^^{^fKso^j^)\^jx). 

If 

" d 






j=l 



then the lifting Z^°°> is defined by the formula (see [4, 5]) 

5 

Z'^°-)=Y,ZW, + 3^^z)- (21) 

It follows from this formula that 



Zi^) = ^ ZW'^ + 3j(^) , (22) 

where 



d sr^ sr^ i d _ s, v-^ v^ „ / , i . „x X 5 



0<|cr|<A: i=l '^ Q<\a\<ki=l '^ 

Let / be the generating function of the contact vector field Z (see formula 
(2)) and Oi = {x,y, z,p,q,u^ ,ul.,Uy,u],,Up,Ug). Then the vector ip{Zf){6i) 
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IS 



{ip^, . . . ,tp^) with 
V'^ = -ulf -ulf:c-u\fy + {-pul-qu\ + u^)fz 

+ {ul + puDfp + {ul + quDfq + fa:x + '^Pfxz + P^fzz 
+ 2v}fxp + (n^ + U^)fxq + 2pU^f^p + p(ti2 + U^)f-^q 
+ {U^ffpp + {U^ + U^Wfpq + V?U^fqq, 

V'^ = -^2/ - n^/^ - ^2/^ + {-pul - qui + n^)/^ 

+ {ul + pu2)/p + {u'i + ^u^)/, + /^y + qf^^ + p/y^ + pqf^^ 

+ ti^/xp + ti''/a;g + « /yp + ^^/yg + {QU^ + PU^)fzp 

+ (gn^ +pu'^)fzq + U^U^fpp + {U^U'^ + {U^f)fpq + U^U^ fqq, 

V'^ = -ti^/ - n^/^ - u^fy + (-pu^ - gn^ + u^)fz 

+ {ul + puDfp + {ul + qul)fq + fxy + qfxz + Pfyz + PQIzz 

+ U^fxp + U'^fxq + v}fyp + U fyq + {QU^ + PU^)fzp 

+ {qU^ +pU^)fzq + U^U^fpp + {U^U^ + {U^f)fpq + u'^U* fqq, 

^^ = -utf - U^fx - U^qfy + {-pUp - qU^q + U^) f . 

/ A . A\ f . / A . A\ fl . fl ^ I 



+ {ui + pn^^fp + {ut + qu^^fq + fyy + 2qjy:, + (^f, 
+ (n2 + u3)/j^p + 2u^fyq + g(u2 + n3)/^p + 2(?u4/. 



zz 
zq 
qq- 



(23) 



3.4. Differential invariants. By T we denote the pseudogroup of all con- 
tact transformations of M. Its action is lifted to J'^vr, /c > 0, as it was 
explained above. 

A function (vector field, differential form, or any other natural geometric 
object on J^tt) is a A; th- order differential invariant of F if for any 99 G 
r the lifted transformation ip^ > preserves this object. In this work these 
differential invariants are called also differential invariants (of order k) of 
Monge-Ampere equations or simply differential invariants (of order k). 

Let £ be a Monge-Ampere equation, S'g the section of vr identified with 
£, and / a differential invariant of order k. Then the the value of I on 8, 
is defined as {jkSg_)*{I) and denoted by /g. If a contact transformation / 
transforms £ to £, then, obviously, f^' transforms /g to /g, for any /cth 
order invariant /. 

Differential invariants that are functions are also called scalar differential 
invariants. By A^ we denote the M-algebra of all scalar differential invariants 
of order < k. By identifying Ak with tt^ f^{Ak) C Ai, V/c < /, one gets a 
sequence of inclusions 

AocAiC.cAkC Ak+i c . . . 

The M-algebra A = IJfc^o ^k is called the algebra of scalar differential in- 
variants of Monge-Ampere equations. 

Remark 3.3. It is worth noticing that a scalar differential invariant / is 
completely determined by its values /g on concrete equations £. This obser- 
vation will be used below. 

Let Z be a contact vector field in M and / a differential invariant of order 
k. Then L^^k) {I) = 0, where L stands for the Lie derivative. This means, in 
particular, that /cth order scalar invariants are first integrals of all contact 
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vector fields lifted to J^tt. Obviously, a scalar differential invariant of order 
k is constant on any orbit of the action of T on J^vr. Such an orbit consists, 
generally, of two components, since contact transformations need not be 
orientation preserving (e.g., the famous Legendre transformation x' = p, 
y' = q, z' = xp + yq — z, p' = x, q' = y is not). In other words, the above- 
mentioned first integrals of Z^ •* are, generally, invariant only with respect 
to the unit component of T and will be called almost invariant. Anyway, 
generic orbits of contact transformations and of contact vector fields have 
the same dimension: 

Proposition 3.4. (1) J^vr is an orbit of the action ofT iff k = 0, 1, 

(2) Codimension of a generic orbit of J'^ir is equal to 2. 

(3) Codimension of a generic orbit of J^ir is equal to 29. 

Proof. Let O^ be a generic point of J^vr and Orbg^. the orbit of the action 
of r on J^n passing through Of^. Then codim Orbg^. = dim j'^ir — dim Orbg^ . 
The dimension of Orbgj, is the dimension of the subspace spanned by all 
vectors X^^>{0k) which can be calculated with the help of computer algebra 
using formulas (22) and (23). D 

An immediate consequence of the above proposition is 

Corollary 3.5. (1) The algebra of scalar differential invariants A2 is 

generated by 2 functionally independent invarints. 
(2) The algebra of scalar differential invariants A^ is generated by 29 
functionally independent invarints. 

Differential invariants constructed below come mainly form natural geo- 
metric constructions without saying that these are invariant with respect 
to the full pseudo-group T. Although not impossible, it is quite challenging 
task to obtain first integrals of Z^ ' analytically even for small k. 

4. Differential invariants on J^tt 

The next step to be done is explicit construction of differential invariants 
that generate A2 as a C°^-closed algebra. 

4.1. Base projectors. Let I) be a distribution on M. Denote by T)^^' the 
distribution generated by all vector fields X and [X, Y], \/ X,Y G D. Setting 
D(o) = V, we define ©('"+1), r = 0, 1, . . ., inductively by the formula I»('"+i) = 
(X)(''))(i). 

Lemma 4.1. For a hyperbolic Monge-Ampere equation £ 

dim(©^)(i) = dim(a3|)(i) = 3. 

Proof Let iv G A^{M) and X,Y e D{M) be such that uj{X) = uj{Y) = 0. 
Then, by applying formula dLo{X, Y) = LxiYJ u;) - Ly(X J u) - [X, Y]Juj, 
one easily finds that 

uj{[X,Y]) = -dLo{X,Y). 

If now Lu = Ui and vector fields X,Y £ D^, i = 1,2, are independent, then 
dUi{X,Y) 7^ due to hyperbolicity of £. So, the above formula shows that 
6''i([X, y]) / 0, i.e., that [X, y] does not belong to the Cartan distribution 
on M. So, [X, Y] is independent on X and Y. D 



MONGE-AMPERE EQUATIONS 13 

Restricting ourselves to the generic case only, we assume from now on 
that 

dim(!D|)(2) = dim(23|)(2) ^ 5 _ ^24) 

Suppose that vector fields Xi, X2 generate the distribution T>^^ and vec- 
tor fields X3, X4 generate the distribution D|. The 3-dimensional generic 
distributions {Xi, X2, [Xi, X2]) and (X3, X4, [X3, X4]) intersect along a one- 
dimensional subdistribution !D| = (Xi,X2, [Xi, X2]) n {X3,X4,[Xs,X4\). 
Hence, equation £ generates a direct sum decomposition 

T(M) = !D|e23|e23|. (25) 

This decomposition generates six projections 

7i:T{M)^Vl, i = 1,2,3, 
Vf : T{M) ^ !D*£ © !D| , j = 1, 2 , 
Te ■.T{M) -^Q = V\®V\. 

These projections may be viewed as vector- valued 1-forms. Namely, let 
X5 be a vector field generating T)\. Consider the co-frame {uj^, . . . ,uj^} on 
M dual to the frame {Xi, . . . , X^}, i.e., LOi{Xj) = 5ij. Then 

•Pi = LO'^ (g) Xi + u^ (g) X2 , 
?2=^^^X3 + UJ^(gX4, 

^3 = w^ X5 , (26) 

Tf^=Tj + J'3, i = l,2, 

^6 = ^1 + ^2. 

These vector-valued differential 1-forms are, obviously, differential invari- 
ants of £ with respect to contact transformations. Moreover, according to 
proposition 3.2, the original equation £ is completely determined by each of 
the projectors ^i, T2- 

4.2. Coordinate-wise description of base projectors. In order to find 
local expressions for the above projectors, consider vector fields Xi, . . . ,^4 
given by (18) and use the notation (20), i.e., 



(27) 



ox oz op oq oy oz op oq 

d d , d n d ^^ d d n d a d 

^3 = TT +P7r + ^ TT + ^ V ' ^4 = TT + ?ir + ^ IT + ^ V • 
ox oz op oq oy oz op oq 

The remaining field X^ is defined by the relation 

X5 = A^Xi + \^X2 + k[Xi,X2] = X'X^ + X^Xi + x[X3, X^] . (28) 

A simple computation shows that 



IMICHAL MARVAN, ALEXANDRE M. VINOGRADOV, AND VALERY A. YUMAGUZHIN 
with 

— 2(ti^ + pu^ + uiUp) — (u^ + u^)u^ + u^Up + u^Up j , 

1 / (29) 

A^ = -^ 7 {u^ + u^)x + p{u^ + u^)z + u^iu"^ + u\ 

w' — u-^ \ 

— 2{uy + qu^ + u^Uq) — (n + u )Up + u Ug + u Ug 

provided that X^ is normahzed by the requirement k = 1. 

Brackets of vector fields Xi, . . . , X^ are described by means of the coeffi- 
cients V"-^: 

5 

[Xj,Xi] = y ^ h)j^Xi . 

4 = 1 

Obviously, 6*^ = -blj. 



4.3. Convenient coordinates on J^'tt. Vector fields Xi, i = 1,...,5 in- 
duce vector fields Xj on the bundle J°°tt, uniquely defined by the condition 
f{Se)*Xi = Xi for ah sections Se- Thus, Xi = Di + pD^ + u^D^ + u'^Dr,, 
etc., where Di denote the total derivatives, see Sect. 3.3.3. 

Differential invariants of hyperbolic Monge- Ampere equations constructed! 
bellow are described in terms of the quantities X^^ . . . Xj^6^-. So, we need to 
know all algebraic relations connecting them, at least for h = 0, 1. To find 
these efficiently it is convenient to use a non-standard local chart in J tt. 

Lemma 4.2. Functions 

^h...ih = 3^n • • • '^ihU^^ k< ■■■ <ih, h<k. (30) 

together with functions x'^,u^ constitute a local chart on J^ir. Moreover, the 
standard jet coordinates on J^vr are rational functions of these coordinates. 

Proof. For k = 2 the assertion is verified directly. For k > 2 one can express 
the standard jet coordinates uj ^ = Di^^^^i^u^ in terms of coordinates (30) 
by making use of the following obvious facts. First, fields Di are linear com- 
binations of fields Xj with coefficients in C°°{J'^it). Second, the coefficients 
6j^J2 are functions on J^vr. Third, XjjXj^/ = — ft^^^^Xj/ + Xj^Xjj/ for every 
function / G C°°(J^7r), k>2. D 

A complete system of relations connecting functions 6^- can be found by 
routine computations and taking into consideration geometric properties of 
fields Xi, . . . ,X5. For instance, 6^2 = ^12 ~ 0, since [Xi,X2] belongs to the 
distribution (I'g)^^) generated by Xi,X2 and X5, etc. The final result is as 
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follows: 



bU = 


0, 


bk = 0, 








b\.= 


0, 


bi2 = 0, 


b'u 


= 1, 




bls = 


-bis, 


bis = -bis, 


bis 


= 0, 




bls = 


-bis, 


°23 — "23- 


bis 


= 0, 




bl,= 


-bh. 


b\, = -bl„ 


b\. 


= 0, 




bl.= 


-bk. 


b\, = -hi,., 


bL 


= 0, 




bl,= 


-b\2. 


bi, = -bh. 


bL 


= -1, 




bl,= 


-bl. 


-bis, 


bk = 


-bl,- 


-b\s. 


bh = 


-b\s 


-bis, 


bl,= 


-bl.- 


-b\„ 


bh = 


-bh 


+ bl5 + h\^, 









(31) 



Henceforth we shall simplify the notation by using Xi for Xj. 

4.4. Curvatures. Using formulas (3), (4) and the direct sum decomposition 
(25), it is easy to compute the curvature forms of projectors ^i, ^2, '^i , 
^2 , -^C) which are 

3?i = w^ A tj^ (g) Xs , 

3?2 = -ur" A cj'' X5 , 

3?i = -{hl^oj^ + hl^uj"^) Auj^(^Xs- {bj^uj^ + bi^LJ^) A w^ X4 , (32) 

Jll = -{bl^io^ + bl^io^) Aio^^Xi- {bl^LO^ + bl^Lo'^) Alo^0X2, 

0i = Jli + Jl2, 

respectively. It is clear that these curvature forms are differential invariants 
of £. 

Frolicher-Nijenhuis brackets of base projectors give new invariant vector- 
valued forms. These, however, turn out to be linear combinations of cur- 
vature forms. More exactly, a direct computation, which is omitted, shows 
that 

[?l,T2l = i(-[yi,yil - [y2,5'2l + 1^3, 3^31), 
[Tl,T3l = K-I^l'^ll + 1^2,^2! - 1^3, 0^31), 
1^2, ^31 = ^(I3'l,Tll - 1^2,^2! - [T3,5'3l) 



and 



1^1, ^il = -2(3^1 + Jii), [T2, ^21 = -2(^1 + 3^2) , 
{93,^3} = -2(011 +JI2). 
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4.5. Scalar invariants on J^tt. The following three invariant 5-fornis with 
values in Ti| = {X^): 

i {Oil J Jii) J {Oil JJli) = Ai Lo'^ A ... A LO^ <S) X5 , 

1{JI\J012)J{'31\J^2) =A2W^ A...Ac^^®X5, (33) 

{■Jll J Oil) J {Oil J 3^2) = A12 w^ A . . . A t^^ X5 , 



with 



Ai = bl^bL - bLbl., A2 = btbL - bf.bi 



(34) 



-"35^45 ~ '^35"455 ^*-2 — "15"25 ~ '^15"255 

A12 = fo?5^35 + bf^bk + b^bh + 4^l5 • 
are proportional. Therefore, the corresponding proportionality factors are 
scalar differential invariants. In particular, such are 

':_-'''"'• (35) 

/2=Ai2/A2. 

Below it will be shown that Ai, A2 are nowhere zero. 

Theorem 4.3. The algebra of scalar differential invariants on J^vr is gen- 
erated by the invariants I^ and I^ . 

Proof. In view of Corollary 3.5, it is sufficient to show that I^ and I^ are 
functionally independent (on J^vr). But this is straightforward from (31). □ 



Coefficients A^-, a = 1, 2, 12, introduced in (34) have a geometrical mean- 

3 

£ 



ing explained below. Fix a generator W = fX^ in 1)% and consider maps 



defined by formulas 

nf(Z2)=Ti([Z2,VF]), nf(Zi) = T2([^i,VF]), 

with Zi G T)\ Z2 G ©^ Since Ti(T'|) = 3'2(I>^) = both Uf and Df are 
C°^(M)-linear. This is seen as well from their local expressions 

uY = fb]^J®Xi, i = l,2, J = 3, 4, 

d^ = fh^uj^ (^ X,, i = 3,4, J = 1,2. 

Consider also 2-forms pf : D* x 23* ^ M, i = 1, 2, defined by 

pf ([/„ Vi)W = -RiiU,, Vi), Ui, V, G Vi. (36) 

Then, obviously, pY = i^/f)'^^ Aa;^, pY = —{l/f)u!^ Auj^, so that both are 
volume forms of T)^ and D^, respectively. Moreover, we have 

(□r)*(pr) = /'Aipf, (37) 

tr(nronr)=tr(nronr) = /2Ai2. 

Proposition 4.4. // £ is generic, then functions Ai, A2 are nowhere zero. 

Proof. By genericity condition (24), DJ^ and D^ are surjective, hence Ai, A2 
are nonzero. D 
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4.5.1. Now consider operators Vf = DJ^oD^ and V^ = D^oDJ^ acting 
on T>^ and 23^, respectively. It follows from (37) that 

A2 - /2Ai2A + /%A2 (38) 

is the characteristic polynomial for each of them. Another peculiarity of the 
situation is that Dj^ send eigenvectors of V2 to that of Vi and similarly 
for DY- 

The discriminant of polynomial (38) is 

/%A2(/i/'-4). 

Its sign coincides, obviously, with the sign of 

This proves that generic hyperbolic Monge-Ampere equations are subdi- 
vided into three subclasses as follows: 

(1) subclass "h": the operator Vj^ has two different real eigenfunctions 

^ /l/2(/l/2_4)>Q^ 

(2) subclass "p" : the operator Vj has a unique real eigenfunction 
^ /i/2(/i/2 - 4) = 0, 

(3) subclass "e" : the operator Vj has no real eigenfunctions 
^ /i/2(/i/2_4)<o. 

4.5.2. Some almost invariants. The previous considerations lead to an al- 
most invariant choice of generator W = fX^ in D|. Namely, define functions 
A^ , i = 1,2, by relations 

(□rr(pr)=Ar/>r, iDYnpY)=AVpr- 

Obviously, AJ^ = /^Aj. This shows that, up to sign, vector fields 

1 

do not depend on the choice of W. In particular, A^ '' = Aj, so that 

1 

By duality, 1-forms 

^i = ^A^\^5, « = 1,2, 
are almost invariant as well. 

It is not difficult to construct further almost invariant forms. For instance, 
the forms 

^^j=%J^j, i = l,2, 
are manifestly almost invariant and have the following local expressions: 

T&ij = J\Aj\ll>^ a cj^ 1?2j = \/\^j\^^ A CJ*^ . 
The products 

pj = (-sign Aj)^ij A '&2J = Aj uj^ A uj'^ A uj^ A uj^, j = 1, 2, (39) 



Wi = , W, i = l,2, 



Wi = ^^X^, i = h2. 
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which are volume forms on the Cartan distribution Iig, ©D|, are, obviously, 
fully invariant. This is a very simple example on how an invariant can be con- 
structed from almost invariants. Forms pj can be described in a manifestly 
invariant way as follows: 

where ( •, • ) stands for convolution. Note that the form 

pi2 ^ Pp. ^ ^^2 cj^ A w^ A w^ A w*^ (40) 

is invariant too. 

Similarly, one can construct many other invariant forms. Some of them 
are : 

{'Jl\j-Jll)=- {bl^LU^ + bl^Uj'^) AUJ^ + {bl^LU^ + bl^Lu'^) A LO^ 

{OIIJOI2) =(fe?5tc^ + 6^5^^) Auj^- {bjr^uj^ + 6^5^^) A uj^ 
Ji\j{ 3^2 J 3^1 > = 2A2 cji A Lj2 A io^ 
JIIj{ JiljJli) = 2Ai uj^ Auj^ Aio^ 
Now it is easy to construct almost invariant volume forms : 

-dj A Pj = |Aj|^/^ u;^ A u;^ A u;^ A LJ^ A uj^, j = 1, 2. 



(41) 



(42) 



5. Differential invariants on J^tt 
Since u!'^{Xi) = const, namely, Sm, we have 

diOk{Xi,Xj) = -io''{[Xi,X,]). 
(see the proof of lemma 4.1). This implies the useful formula 



dio" 






UJ 



Alo^ . 



(43) 



5.1. The complete parallelism. First, note that invariant differential 1- 
forms dl^ and dl^ live on J^tt. This leads us immediately to another set of 
invariant differential 1-forms on J^tt: 

Q} = ?i Jd/^ = Xl(/l)c^l+X2(/^)a;^ 
f]2 = Ti J (i/2 = Xi (/2) UJ^ + X2{I'^) u? , 
f]3 = T2Jd/i = X3(/1)cj3 + X4(/^)l<j\ (44) 

f]^ = T2 J d/2 = X3(/2) U? + X4(/2) Lj4 , 

VL\ = '^^AdI^ =X5(/1)CU^ fi^ = ?3Jd/2 = X5(/2)t^5_ 

Supposing that £ is a generic equation, we henceforth assume that 

X5(/i)/0, X5(/2)/0, (45) 

and 



Ai 



Xi(/i) X2(/i) 

Xi(/2) X2(/2) 



/O, A2 



X3(/l) X4(/l) 
X3(/2) X4(/2) 



/O. 



(46) 



This means that two sets of forms {^2^, . . . , fi'^, i7f } and \S^ ■, ■ ■ ■ , 0,^, ^2} 
are invariant coframes on M (we omit the subscript £ according to Remark 
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3.3). Each of these coframes determines an invariant complete parallehsm 
on M. 

The frames {Yi, . . ., Y4, Y^} and {Yi, . . . , I4, Y^}, dual to the above con- 
structed coframes, are, obviously, invariant. An explicit description of them 
is : 

Y, = ^{X2{I^)X^ - X^(I^)X2) , 

Y2 = -^{-X2{I^)Xi + Xi{I^)X2) , 

^3 = -^(X4(/2)X3 - X?XI^)X^) , (47) 

Y^ = -L(-X4(/i)X3 + X^{I^)X^) , 

^2 

5.2. More scalar invariants on J^vr. Among numerous invariants con- 
structed previously there are functions, (vector-valued) differential forms, and 
vector fields. Further invariants can by obtained just by applying various op- 
erations of tensor algebra, Frolicher-Nijenhuis brackets, etc, to these objects. 
Moreover, components of an invariant object with respect to an invariant 
basis are scalar differential invariants as well as its proper differential invari- 
ants. These simple general tricks are rather efficient and were already used 
in constructing differential invariants on J^vr. As for J^tt we shall proceed 
along these lines as well. 

The invariant 1-forms ilf and Ol are proportional. So, the proportionality 
factor 



,3 _ MI'] 



(48) 



Mn 

is a scalar differential invariant on J^vr. 
Consider now invariant 2-forms on J^vr: 

Oil J di^ = i^n^ An'^, 
oiijdi^ = fn^ aq"^, 
oi2Jdi^ = i^n^ aq*, 

oi\jdi^ = i^^n^ Afil + i^^n^ Anl, 
oil J di^ = i'^^' Anl + i^^n^ A nf , 
oil J di^ = i'^^^ Anl + i^^n^ A nf , 
oil J di^ = i'^^^ Anl + i^^n^ A nf . 

Their components /^, . . . , l"^"^ with respect to the base ^^, . . . , i7^ are fur- 
ther scalar differential invariants on J^vr. The simplest among them are 
/6 = Ai/X5(/i) and /8 = A2/X5(/i). 

In the same manner one easily find numerous non-scalar differential in- 
variants on J^TT. For instance, such are 3-forms [yj,3?j] or [yjjDiJ], 4-forms 
Pi, {01]J011)1 5-forms [?„ Ol]\ J pk,Ol]l etc. 
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5.3. Better manageable invariants. From the above said one can see 
that there are sufficient resources for constructing differential invariants and 
the main problem becomes to select functionally independent ones in the 
simplest possible way. From technical point of view this forces us to look for 
manageable invariants, for instance, those that have local expression as sim- 
ple as possible. In the considered context a help comes from almost invariant 
objects as it is illustrated below. 

In view of (39), (40) and (43), for a = 1,2,12 we have the invariant 
5-forms 

dpa = d(A„ Lu'^ Auj'^ Auj^ A w") 

= {X5{A„)+A„B)lu'^ Alo^ Auj^ Auj'^AuJ^, 

where B = bl^ + b^^ + b^^ + bj^ = 2{b\^ + ^2^) = 2(63^ + bj^) according to 
(31). 

By comparing these 5-forms with (42) we obtain almost scalar invariants 

• X5{A^)+AaB 1 o in • 1 o fr.1\ 

^" |A-|3/2 ' 0- = 1,2, 12, J = 1,2, (51) 

on J^vr which are better manageable in comparison to those constructed in 
the previous subsection. The squares (li)"^ are, obviously, full scalar invari- 
ants. Apart from the obvious relation (/^//2)2 = (/^//^^s ^j^gy g^j.g function- 
ally independent. Some of the earlier constructed invariants can be expressed 
in terms of /^'s, e.g., 

Mil = (I'u-Hm' . ^ . 2 
x^iP) [ii^-iiP)p' ^ ' ■ 

6. The equivalence problem 

So far we obtained two independent second-order scalar invariants 1^,^ 
(see (35)) and a number of third-order invariants. Put (see (51)) 

I' = m?, I' = {ii)\ I' = (ihf, 

The following statement can be checked by a direct computer-supported 
calculation in coordinates (30): 

Theorem 6.1. Invariants I^ , P, I^, P, and P are functionally indepen- 
dent on J"3(vr). 

Of course, this choice of basic scalar invariants is not unique. For instance, 
invariants 1^,1^, I"^,!^, I^ (see (48), (49)) are functionally independent as 
well. However, this and other reasonable choices are "less manageable" with 
respect to those made in the above theorem. Unfortunately, this fact is not 
clearly seen from the above exposition, since we were forced to skip technical 
details of computations. 

According to "the principle of n invariants" [22], any quintuple of func- 
tionally independent scalar invariants gives a solution of the equivalence 
problem for generic hyperbolic Monge-Ampere equations. Theorem 6.1 guar- 
antees existence of a such one, namely, P, . . . , P. 
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More exactly, let £ be a generic hyperbolic Monge-Ampere equation con- 
sidered as a section of the bundle tt. Since invariants I^,. . . ,1^ are func- 
tionally independent their values /g, . . . , /| on £ form a (local) chart in M. 
In terms of these coordinates, the 1-forms f^i, . . . , f^s, defining an absolute 
parallelism on M, are described in terms of functions QUl)^, ■ ■ ■ , /|) coming 
from the decomposition 

5 

n, = Y,^^{ih---,4)dii, i = i,...,5. 

Theorem 6.2. The (local) equivalence class of a generic equation £ with 
respect to contact transformations is uniquely determined by the family of 
functions QUlj^, . . . , /|), i = 1, . . . , 5 . 

Proof. Let £ be another generic Monge-Ampere equation such that there 
exists a contact transformation transforming it to £. Then, obviously, the 
functions ^U/g, . . . , /|) and ^){lh ■■■,!?■) coincide for all i and j. 

Let £, £ be Monge-Ampere equations such that for all i and j the func- 
tions ^]{l\, ■■■ -,!%) and ^]{ll-, ■■ ■ ,ll) coincide. Let Ig = {ij-, . . . , l|) and 
/g = (/i ...,/5) be invariant coordinate systems in M for £ and £ re- 
spectively. Then /-T o /g is a locally defined diffeomorphism M -^ M. 
This diffeomorphism is a contact transformation because it transforms fij = 
E •=! ^iilh ■■■, 4)dli to ZU m^l ■■■, Il)dll = f^^ i = 1, . . . , 5, and, 
in particular, the contact form Vl^ to the contact form ils. By obvious reasons 
it also transforms the pair of distributions (!D£,D|) to the pair (Di,!)?) 

and hence £ to £. D 

7. Examples 

Examples discussed in this section aim to illustrate the character and com- 
plexity of problems related with actual computations and use of differential 
invariants. Henceforth invariants I* are denoted by /j. 

Example 7.1. Consider the equation 

'^\^xx^yy ^xy) ' V ^xx ^^V^xy ~r X Zyy -\- X y Z — U. 

The first two invariants are /i = znj^/d, I2 = zn^/d, where 
n± = 2(z + 37/ T 2)x^zl -{z + I2x^y^)xyz^zy + 2{z + 3x^ ± 2)y^z^ 
+ {z^ + Sx^y^z + Ay'^z ±4.z± IQx^y^ ± 16/ - I2)xz^ 
+ {z^ + Ax^z + Sx^y'^z =F 4z =F 16x'' =F lOx^y^ - I2)yzy 
+ 2z^ + mx'^y'^z^ + Qy^z^ - Ax^y^z^ + Qx'^z^ 
- 8z =F 16x^z ± 16/z + 8a;^ + l&x^y"^ + S/, 
d = A{z^ + Zy^z + A)x^zl - 2{z^ + Ux'^y^z - 16)xyz^Zy 

+ A{z^ + ?>x'^z -F A)y^zl -F 2[z^ -F Sx^y^z -^ Vz -^ 2^)xzz^ 
-F 2(z^ -F Ax'^z -F Sx^y^z + 20)yzzj, -F 4(18x''/z^ -F z^ 
-F Sx^z^ - 2x^y^z^ -F S/z^ -^ 12z -F 4x^ -F Sx^y^ + Ay'^)z. 
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The invariants 1^,14, I^ are large fractions whose non-reducible numerators 
are polynomials of order three in Zx, Zy, five in z, and six in x, y. Invariants 
Is, s > 5, are even more cumbersome. 

Computation shows that the jacobian d{Ii, I2, I3, I4, l5)/d{x,y, z, Zx, Zy) 
is nonzero, hence the first five invariants are functionally independent and 
can be chosen to be local coordinates on J^{t). Although an explicit in- 
version is rather hopeless, one can still find algorithmically the relations 
connecting principal invariants Ii, . . . ,1^ and higher I^ at least in principle. 
This kind of procedure is outlined in Example 7.2 below. 

Example 7.2. Put ^ = Zx + Zy + e and consider the family of equations 

(iZxZy + C'^){ZxxZyy - Zly) + AC'^{ZyZxx + ZxZyy + C^) = 0. (52) 

depending on parameter e. Assuming that e 7^ 0, we have 

_ {zx + Zyf + "Jejzx + Zy) + 4e^ 

hezx + ezy + 4e2 
_ (zx + Zyf + 3e(z^ + Zy) + 4e^ 

ezx + bezy + 46^ 

^^ " ^ ei/2(5z^ + z^ + 4e)3/2 ' 

3/2 ~^x + ^ZxZy + 7Zy + Sezi: + 33ezy + 21e^ 
^^ " ^ ei/2(5z^ + ^j^ + 4e)3/2 ' 

^ 5/2 (^^ + %)^ + Mzx + %)^ + 17e^(2^ + Zy) + 216^ 
^ e3/2(5z^ + z^ + 4e)3/2 

All invariants are independent of x, y, z, reflecting the fact that x ^^ x + ti, 
y *-^ y + t2, z \-^ z + t3 are symmetries of equation (52). One easily checks 
that /i,/2 are functionally independent, but it is still not straightforward to 
express Zx, Zy in terms of /i, I2 explicitly. 

To establish the dependence of I^, s > 3, on hjh, we observe that for 
every s there exists a polynomial Ps{zx,Zy,Is) such that Is is a solution 
of the equation Pg = 0. Then what we need is eliminating Zx,Zy from the 
system 

^izx + Zy)"^ + 2>e{zx + Zy) + 4e2 
-'I ~ ^ ;: ; ~~ro = '-'' 



f>ezx + eZy + Ae^ 

, {zx + Zy)'^ + 3e{zx + Zy) + 4e^ 

eza: + bezy + 4e2 



0, 



^syZxt Zy, Is) — U. 



To this end, it suffices to compute the Grobner basis of the last system with 
respect to an "elimination ordering" of monomials. With the help of the 
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Groebner package of Maple 10 the following quadratic equation for I^, 

= 4096 /|/| - /|(729 /f /| - 1971 /f /| + 20493 llll 

+ 3563 if /2 - bllUllll + 183915 /i/| 

+ 3951 if - 52723/1^/2 + 45517 /i/| + 102191 ll)h 

+ (27/f /| - 81 /f /| - 32 Ifh + 426 Ipl - 1206 /i/| 

- 44 /f + 270 /f/2 - 800 /i/| - llUllf 

can be found rather quickly as well as similar quadratic equations for I^^I^. 
The assumptions of Sect. 5.1 are satisfied as well. In particular, Ai, A2 are 
nonzero since 

, {zx + Zy){^Zx + ^Zy + Se){z,j, + Zy + c)^ 



Ai = A2 = -128- 



X 



e^{zx + 5zy + 4e)2(5z^ + Zy + 4e)2 
z"^ + 2zxZy + Zy + 3ezx + 3ezy + 4e^ 
z^ + 622:^2/ + ^y + 2eza: + 2ezj^ + e^ 



This enables us to compute the higher invariants. For instance, Iq is solution 
of the quadratic equation 

= -16lf{27lfl2 - 27/f/| + 22/f - 56/f/2 - 2/f/| 

+ 8/1^/2 - 42/f + 50/1 /| + 28/1^ + 56/1/2 + 28/|)/| 

+ /i(/i/2 - /l - /2)(9/i/2 + 7/1 + 7/2)(3/f/2 - 3/f/| 

- 26/f - 34/f /2 - 8/1 /| + 18/1^ + 36/1/2 + 18/|)/6 

+ (/i + hfihh -h- h){^hh + 7/1 + 7/2)(/i/2 - 2/1 - 2/2)2. 

Although every invariant computed so far depends on e, its expression in 
terms of /i,/2 does not. This suggests the idea that the parameter e is 
removable. And indeed, after substitution z ^^ ez equation (52) becomes 
equivalent to itself with e = 1. Thus, the family of equations (52) consists 
of a continuum of generic members with e 7^ 0, which are all mutually 
equivalent, and a single non-generic member with e = (in which case 
Ai = A2 = 0). 

Example 7.3. Consider the family of equations 

"gyZxxZyy Zxy) + y Zxx ^xyzxy + X Zyy + ex y — u, 

depending on a real parameter e 7^ 4. Then the first five invariants are 
constants 

J _j . .e + 12 
ii — H - 



h = 
3200- 



e-4' 
800 



e-4' 
[e + 12)2 



(e - 4)3 ' 

while the higher invariants Is are undefined. 

The equation belongs to the subclass "h", or "p", or "e" (see 4.5.1) if 
e > —4 or e = —4 or e < —4, respectively. 



2michal marvan, alexandre m. vinogradov, and valery a. yumaguzhin 

Acknowledgements 
M. Marvan acknowledges the support from GACR under grant 201/04/0538. 

References 

[1] D.V. Alekseevskiy, A.M. Vinogradov and V.V. Lychagin, Basic ideas and concepts of 

differential geometry, in: Geometry, I Encyclopaedia Math. Sci. 28, Springer, Berlin, 

1991, 1-264. 

[2] A. Frolicher and A. Nijenhuis, Theory of vector valued differential forms. Part I: 

Derivations in the graded ring of differential forms, Indag. Math. 18 (1956) 338-359. 

[3] P. Hartman and A. Wintner, On hyperbolic partial differential equations, American 

Journal of Mathematics 74 (1952) 834-864. 
[4] I.S. Krasil'shchik, V.V. Lychagin and A.M. Vinogradov, Geometry of Jet Spaces and 

Nonlinear Partial Differential Equations, Gordon and Breach, New York, 1986. 
[5] I.S. Krasil'shchik and A.M. Vinogradov, Editors, Symmetries and Conservation Laws 
for Differential Equations of Mathematical Physics, Translations of Mathematical 
Monographs. Vol.182, Providence RI: American Mathematical Society, 1999. 
[6] B.S. Kruglikov, Some classificational problems in four-dimensional geometry: distri- 
butions, almost complex structures and the generalized Monge-Ampere equations. 
Math. Sbornik 189 (1998) (11) 61-74 (in Russian); English translation in Sb. Math. 
186 (1998) (11-12) 1643-1656; e-print: http://xxx.lanl.gov/abs/dg-ga/9611005. 
[7] B.S. Kruglikov, Symplectic and contact Lie algebras with application to the Monge- 
Ampere equations Trudy Mat. Inst. Steklova 221 (1998) 232 — 246 (in Russian); 
English translation in Proc. Steklov Math. Inst. 221 (1998) (2) 221-235; e-print: 
http://xxx.lanl.gov/abs/dg-ga/9709004 
[8] B.S. Kruglikov, Classification of Monge-Ampere equations with two variables, in: Ge- 
ometry and Topology of Caustics - CAUSTICS '98 (Warsaw), Banach Center Publi- 
cations 50 (Polish Acad. Sci., Warsaw, 1999) 179-194. 
[9] A. Kushner, Monge-Ampere equations and e-structures, Dokl. Akad. Nauk361 (1998) 
(5) 595-596. 

[10] H. Lewy, Uber das Anfangswertproblem bei einer hyperbolischen nichtlinearen par- 
tiellen Differentialgleichung zweiter Ordnung mit zwei unabhangigen Veranderlichen, 
Math. Annalen 98 (1928) 179-191. 

[11] V.V. Lychagin, Contact geometry and non-linear second order differential equations, 
Russian Math. Surveys 34 (1979) 149-180. 

[12] V.V. Lychagin, Lectures on Geometry of Differential Equations, Universita "La 
Sapienza," Roma, 1992, 133 p. 

[13] V.V. Lychagin, V.N. Rubtsov and I.V. Chekalov, A classification of Monge-Ampere 
equations, Ann. Sc. Ecole Norm. Sup. (4) 26 (1993), 281-308. 

[14] M. Marvan, A.M. Vinogradov and V.A. Yumaguzhin, Differential invariants of generic 
hyperbolic Monge-Ampere equations Russian Acad. Sci. Dokl. Math. 405 (2005) 299- 
301 (in Russian) English translation in: Doklady Mathematics 72 (2005) 883-885. 

[15] M. Matsuda, Two methods of integrating Monge-Ampere's equations. Trans. Amer. 
Math. Soc. 150 (1970) 327-343. 

[16] M. Matsuda, Two methods of integrating Monge-Ampere's equations. II, Trans. 
Amer. Math. Soc. 166 (1972) 371-386. 

[17] T. Morimoto, Monge-Ampere equations viewed from contact geometry, in: Symplectic 
Singularities and Geometry of Gauge Fields (Warsaw, 1995), 105-121, Banach Center 
Publ., 39, Polish Acad. Sci., Warsaw, 1997. 

[18] O.P. Tchij, Contact geometry of hyperbolic Monge-Ampere eqquations, Lobachevskii 
Journal of Mathematics 4 (1999) 109-162. 

[19] D.V. Tunitskiy, Monge-Ampere equations and functors of characteristic connection, 
Izv. RAN, Ser. Math. 65 (6) (2001) 173-222. 

[20] D.V. Tunitskiy, On the global solvability of hyperbolic Monge-Ampere equations, 
Izv. Ross. Akad. Nauk Ser. Mat. 61 (1997), No. 5, 177-224 (in Russian); translation 
in Izv. Math 61 (1997), No. 5, 1069-1111. 



MONGE-AMPERE EQUATIONS 25 

[21] A.M. Vinogradov, Geometry of nonlinear differential equations, manuscript, DiSiety 

Institute, 1987. 
[22] A.M. Vinogradov, Scalar differential invariants, diffieties and characteristic classes, 

in: Mechanics, Analysis and Geometry; 200 Years after Lagrange, ed. M. Francaviglia 

(North- Holland), pp. 379-414, 1991. 
[23] A.M. Vinogradov and V.A. Yumaguzhin, Differential invariants of webs on 2- 

dimonsional manifolds, Mat. Zametki 48 (1990), No. 1, 46-68 (in Russian). 

Mathematical Institute of the Silesian University in Opava, Na Rybni'cku 1, 
746 01 Opava, Czech Republic 

E-mail address: Michal . MarvanOmath . slu . cz 

University of Salerno, Dipartimento di Matematica ed Informatica, Fa- 
colta' di Scienze Matematiche, Fisiche e Naturali, Universita' degli Studi di 
Salerno, Via Ponte del Melillo, 84084 Fisciano (SA), Italy, and Istituto Ital- 
lANO di Fisica Nucleare, Italy 

E-mail address: vinograd@unisa.it 

Program Systems Institute of RAS, 152020, Pereslavl'-Zalesskiy, m. Botik, 
Russia 

E-mail address: yuma@diffiety.botik.ru 



